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ABSTRACT 


To  predict  motions  of  SWATH  ships  in  following  seas, 
especially  when  the  encounter  frequency  is  small,  unified 
slender-body  theory  has  been  applied.  The  longitudinal 
interaction  term  is  computed  and  added  to  the  results  of 
the  strip  theory  as  a  correction  term.  The  added  mass 
coefficients  are  computed  to  be  much  larger  than  those 
from  strip  theory,  while  damping  coefficients  are 
slightly  less  than  two-dimensional  results.  There  is  im¬ 
provement  for  heave  motion  for  a  limited  range  of  fre¬ 
quencies,  but  the  pitch  motion  is  generally  poor.  To 
improve  these  results  further,  the  two-dimensional  potential 
should  be  solved  by  the  multipole  expansion  method  rather 
than  by  the  Frank  close-fit  method. 


ADMINISTRATIVE  INFORMATION 

This  study  was  performed  under  the  General  Hydromechanics  Research  Program  and 
was  authorized  by  the  Naval  Sea  Systems  Command,  Hull  Research  and  Technology  Office. 
Funding  was  provided  under  Program  Element  61153N,  Task  Area  SR  0230101,  and  Work 
Unit  1572-031. 


INTRODUCTION 

An  analytical  method  to  predict  the  motions  of  SWATH  ships  has  been  developed 

1* 

by  Lee.  The  numerical  results  computed  by  this  method  provide  good  correlation 

with  model  test  results  for  moderate  speed  ranges.  The  present  author  has  improved 

the  prediction  of  heave  and  pitch  motions  in  head  seas  by  adding  surge  effect  to 

2 

the  pitch  exciting  moment  and  by  correcting  the  viscous  damping  terms. 

However,  in  following  seas  correlation  between  the  analytical  method  based  on 
strip  theory  and  model  test  results  is  not  satisfactory.  In  particular,  when  the 
SWATH  ship  moves  almost  as  fa:  t  as  the  wave  celerity,  the  encounter  frequency  be¬ 
comes  very  small  and  application  of  strip  theory  is  not  valid.  The  fundamental 
assumption  of  strip  theory  is  that  the  frequency  is  far  larger  than  the  product  of 
the  longitudinal  gradient  of  the  body  surface  and  the  forward  speed. 

3 

Newman  has  recently  applied  slender-body  theory  to  the  problem  of  ship  motions. 
He  has  developed  the  unified  slender-body  theory  which  is  valid  for  all  frequencies. 


*A  complete  listing  of  references  is  given  on  page  75, 
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In  this  theory  the  longitudinal  interaction  terra  is  computed  by  matching  the  inner 

and  outer  solutions.  This  term  and  the  results  of  the  strip  theory  encompass  the 

solution  of  the  unified  slender-body  theory. 

The  numerical  results  for  added  mass  are  larger  than  those  from  strip  theory 

alone  and  the  damping  coefficients  are  generally  smaller  than  the  results  from  strip 

theory.  The  heave  and  pitch  motions  generally  compare  well  with  experimental  result 

when  the  encounter  frequency  is  small.  For  high  encounter  frequencies,  the  results 

become  extremely  large,  especially  for  pitch  motion. 

The  large  discrepancies  at  high  frequencies  might  be  explained  by  the  fact  that 

in  solving  the  two-dimensional  potential  (strip  theory),  we  have  applied  the  Frank 

3 

close-fit  method,  while  Newman  applied  the  method  of  multipole  expansion.  In  multi 
pole  expansion,  there  is  a  clear  separation  of  the  source  strength  at  the  origin 
from  that  of  the  wave-free  potential  while  in  the  Frank  close-fit  method  there  is  no 
separation.  These  different  approaches  to  the  solution  of  the  two-dimension  poten¬ 
tial  function  is  believed  to  cause  these  large  discrepancies. 

BOUNDARY -VALUE  PROBLEM 

We  define  two  coordinate  systems:  the  first,  0  X  Y  Z  is  fixed  in  space  and 

oooo  1 

tlie  second,  oxyz,  is  fixed  with  respect  to  the  ship  wbicli  moves  with  forward  speed, 

U,  in  the  positive  O^X^-axis.  The  oz-axis  is  directed  vertically  upward  and  the 
ox-axis  is  positive  in  the  direction  of  the  ship's  forward  velocity;  see  Figure 
1.  The  oxy-plane  is  the  plane  of  the  undisturbed  free  surface.  These  two  co¬ 
ordinate  systems  coincide  when  the  ship  is  at  rest. 


Figure  l  -  Coordinate  System 


Let  the  surface  of  the  ship  be  specified  mathematically  by  the  equation 


y  =  h(x, z , t)  (1) 

and  let  the  free  surface  be  given  by 

z  =  C(x,y,t)  (2) 

Then,  the  fluid  motion  can  be  expressed  by  the  velocity  potential  <}>(x,y,z,t) 

4>(x,y,z,t)  =  (J)  (x,y,z)  +  <j)2(x,y,z)e  1UJt  (3) 


where  <j>^  =  -Ux  +  (J>o(x,y,z),  the  potential  due  to  steady  forward  motion,  and  <t>^  =  un¬ 
steady  potential  due  to  oscillation. 

Equation  (3)  must  satisfy  the  following  conditions: 

1.  The  Laplace  equation  in  the  fluid  domain: 


4)  +4’  +4'  =0 

xx  yy  zz 


(A) 


2.  The  dynamic  free-surface  condition: 

i'  +  g c  +  7  (42+4>2+,J)2)  =  4  U2  for  z  =  r, ( x ,  y ,  t )  (3) 

l  z  x  y  2  z 


3.  The  kinematic  free-surface  condition: 


4  s  +  4>  C  -  4>  +  t  =  0  for  z  =  <;(x,y,t)  (6) 

x  x  y  y  z  t 

4.  The  kinematic  body  condition: 

4'h  -  (J)  +  4>h  +  h  =  0  for  y  =  h(x,z,t)  (7) 

xx  y  zz  t 
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5.  The  radiation  condition;  that  is,  the  energy  flux  of  waves  associated  with 
the  disturbance  of  the  ship  is  directed  away  from  the  ship  at  infinity, 
liquations  (5),  (6),  and  (7)  are  exact  boundary  conditions.  By  substitution  of 
liquation  (3)  into  Equations  (5)  and  (6),  and  hy  linearization,  the  free-surface 
condition  is  given  by 


U,  ,  .  .  I  .  oj  2i(nLT  .  U  ,  l  -icot  n 

—  +  <t>  +  {  - -  •!>,,  +  —  l|>,,  /  e  =0 

g  oxx  oz  |  2z  g  2  g  2x  g  2xx  | 


i'he  body  condition  can  be  expressed  as  follows: 


(V<{:  +V*„ +U  )  •  n  =  U„ 
o  2  o  2 


where  II  =  (-11,0,0),  steady  forward  velocity  vector 
•> 

U9  =  velocity  vector  due  to  oscillation 

n  =  unit  normal  vector  directed  into  the  fluid  domain 
The  solution  of  Equation  (4)  with  two  boundary  conditions,  Equations  (8)  and  (9), 
is  nonlinear  between  and  %*>__.  The  usual  practice  for  solving  this  problem  is  to 
separate  t>°  and  and  evaluate  these  two  potential  functions  independently.  By 

substitution  of  Equation  (3)  into  Equation  (4)  and  by  separation  of  and  <p ^  in 
Equations  (8)  and  (9),  we  have  two  sets  of  partial  differential  equations. 


|  —  0 
oxx  oyy  ozz 


—  ;  +  •;  =  o 

g  oxx  oz 


(V;  +11)  •  n  =  0 
o  o 


2xx  '2yy  ' 2zz 
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,  w2  ,  2i toU  .  u2  ,  . 

<t>  -  -  lb  +  -  d>_  H - d>_  =0 

2z  g  2  g  y2x  g  v  2xx 


(14) 


W2-U2)  •  n  =  0 


(15) 


Equations  (10),  (11),  and  (12)  comprise  the  potential  function  due  to  the  steady 
forward  motion  and  Equations  (13),  (14),  and  (15)  comprise  the  potential  functions 
due  to  oscillation  with  steady  forward  motion. 

VELOCITY  POTENTIAL  OF  STEADY  FORWARD  MOTION 

The  solutions  of  Equations  (10),  (11),  and  (12)  lead  to  Michell's  integral  in 

the  thin-strip  wave  resistance  problem.  In  the  context  of  slender-body  theory, 

however,  this  three-dimensional  problem  becomes  a  two-dimensional  problem  which  has 

4-6 

been  discussed  in  depth  by  Tuck.  The  solution  of  two-dimensional  potential  is 

given  by 


log 


r  +  |  log 


r^  -  Ux 


(16) 


where  a  =  ^  S';  the  strength  of  source 

S  =4  Trr2;  area  of  immersed  cross  section 
2  o 

—2  2  2 

r  =  (y-nr  +  (Z -rX 

—2  2  2 

r“  =  (y-n)  +  (z+r,) 

Here,  (y,z)  is  the  point  where  the  potential  is  solved  and  (n,0  is  the  source  point. 

UNSTEADY  POTETIAL  DUE  TO  OSCILLATION 

The  solution  of  Equations  (13),  (14),  and  (15)  is  a  three-dimensional  problem. 
This  can  be  solved  by  distribution  of  the  three-dimensional  source  strength  over 
the  surface  of  the  ship  and  by  numerical  solution  of  the  integral  equation.  How¬ 
ever,  this  method  requires  considerable  computer  time  in  numerical  integration.  To 

3 

avoid  such  a  lengthy  process,  the  slender-body  theory  has  been  introduced  by  Newman 
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and  Newman  and  Tuck7  to  solve  this  problem.  Because  the  unsteady  motions  are  assumed 
small,  the  potential  <f>2  in  Equation  (3)  can  be  decomposed  linearly 

j  =  6 

*2  ‘  *>„  +  *7  *£  5j  *1  ll7> 

j  =  l 

where  =  incoming  wave  potential 

o 

=  diffraction  potential 

=  velocity  potential  due  to  motions  of  the  ship  with  unit  amplitude  in 
j  each  of  six  degrees  of  freedom 
r, .  =  amplitude  of  motion  in  each  of  six  degrees  of  freedom 

The  diffraction  potential,  *>7,  satisfies  the  following  condition 


f-  (  *  +  *7>  =  0 

o  / 


for  y  =  h(x,z,t)  (18) 


and  the  v 


elocitv  potential,  * .  (j-1,2, . . .6),  satisfies  the  body  boundary  conditions 


'fi  =  -  ion  .  +  m  . 
j  n  j  J 


Uvre,  the  compon 


ents  of  the  unit  vector  are  defined  as 


(nf ,n9 ,n2)  n 


(n4’n5’n6)  =  (xXn) 


i.  are  defined  by  Ogilvie  and  Tuck  as 


(m^m^m^)  =  m 


id  =  -  (  n  •  V )  Vih 


6 


(23) 


(m^,m5,m6)  =  -(n*V)  (x^V^) 


where  x  is  a  position  vector. 

In  the  slender-body  theory,  the  length  of  the  ship  is  assumed  far  larger  than 
the  beam  and  draft.  With  this  assumption,  the  components  in  Equations  (21),  (22), 
and  (23)  reduce  to 


hh 


+n 


-xn3,  xn2) 

(24) 

(25) 

(26) 

By  substitution  of  Equation  (17)  into  Equation  (14),  the  free-surface  condition  for 
the  potential,  ^  (j=0,l, . . . 7) ,  is  given  by 


g 


jz 


-  0) 


$ .  +  2 iwU  y .  + 

J  jx 


U2  <p. 
jxx 


=  0 


for  z  =  0 


(27) 


With  the  assumption  of 


S 


jz 


-  OJ 


V.  = 
J 


for  z  =  0  (28) 


Equation  (27)  will  be  applied  in  the  outer  region  where  the  three-dimensional  so¬ 
lution  is  expected  and  Equation  (28)  will  be  applied  in  the  inner  region.  The  two- 
dimensional  Green  function  which  satisfies  Equation  (28)  is  given  by  Wehausen  and 

9 

Laitone  for  a  single-hull  body  as 


G2D(y’z)  =  7  PV  J  e-~K— -°k~ky  dk  "  ^  COS  Ky  (29) 

o 
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Equation  (29)  is  the  potential  function  of  a  unit  source  at  the  origin  where  K  = 

2 

.“/g.  With  a  change  in  the  integral  path,  the  Green  function  can  be  expressed  as 


G2D(y’z)  =  ¥  { 


kz 

e  cos 
K  -  k 


dk 


(30) 


The  contour  of  the  integral  path  is  around  the  pole  from  beneath.  With  a  change  in 
the  variables  Equation  (30)  can  be  written  as 


G 


2D 


-  Re  {eKz+lKyE1  (Kz+iKy)} 

7T  1 


.  Kz+iK|y| 
le  1  J  1 


(31) 


The  derivation  of  Equation  (31)  is  given  in  .Appendix  A.  The  exponential  integral  is 
K.(u)  and  is  defined  as 


E1(u)  = 


(32) 


The  asymptotic  properties  of  the  two-dimensional  Green  function  can  be  obtained  from 
the  corresponding  approximation  of  the  exponential  integral.  For  small  values  of  Kr, 
Equation  (31)  can  be  expressed  as  (see  Appendix  B) 

°2D  =  [Y+9-n(Kr)-i7T]  (33) 

where  ,  =  0.577...  which  is  Euler's  constant  and  (r,G;  are  polar  coordinates  such 
th.it  v  =  r  sin  9  and  z  =  -r  cos  0.  For  large  values  of  K  ]  y  |  the  asymptotic  approxi¬ 
mation  of  Equation  (31)  represents  the  outgoing  two-dimensional  plane  waves  in  the 
form 


G 


2D 


leKz+iK|y| 


for  K | y |  >>  1  (34) 
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The  three-dimensional  Green  function  which  satisfies  Equation  (27)  is  given  by 

g 

Wehausen  and  Laitone  as 


oo  n 

G(x,y,z)  =  rl  f  du  f  af.-e.2U.exj.,0:ixu  cos  9-iyu.sin  9_) 

4tt  J  J  gu  -  (art-Uu  cos  0) 

f\  -TT  c' 


de 


(35) 


Equation  (35)  is  the  potential  function  of  a  unit  source  which  is  located  at  the 
origin  and  moves  with  constant  velocity  l).  Equation  (35)  is  multiplied  by  a  constant 
value,  -1/4tt,  for  convenience.  With  the  change  of  variables,  u  cos  0  =  k,  u  sin  0=4 
and  ududG  =  dkd4,  C(x,y,z)  can  be  rewritten  as 


G  (x 


,y,Z)  =  -  j*  e'lkx  dk  -1-  J 


_s* 


2  2 
z(k  +4  ) 


1/2 


-iy4 


g(k2+42) 


1/2 


d4 


-  (co+kU)' 


(36) 


If  we  define  the  Fourier  transform 


f*(k) 


o-> 

■  |  f(;l 


\  ikx  j 

)  e  dx 


(37) 


and  the  inverse  Fourier  transform 


f  (x) 


OO 

=  k  J 


f*(k)  e  lkx  dk 


(38) 


then  the  Fourier  transform  of  the  three-dimensional  Green  function  is  given  by 


O 

G*(y,z;k)  =  -  ~  J 


2  2  1/2 

z(k"+«/)  -iy 4 

ge _ e 

1/2 


d4 


g(k2+e2)  -  (w+kU) 2 


(39) 
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The  value  of  G*  for  k  =  0  reduces  to  the  two-dimensional  source  potential.  Equation 
(30),  as 


G*(y,z;0) 


1_ 

2tt 


CO 


ez|S.|e-iyH 

U!  -  K 


d  £ 


(40) 


An  approximation  similar  to  Equations  (33)  and  (34)  can  be  derived  for  the  transform 
G*(y,z;k).  An  asymptotic  expansion  of  Equation  (39)  for  Kr  <<  1  is  derived  by 
Ursell^  in  the  special  case  U  =  0. 

For  the  case  U  >  0,  the  asymptotic  expansion  is  given  by  (see  Appendix  C) 


G*(k,x)  =  G2D  -  j  f*(k,K,K) 


(41) 


where  G^D  is  Equation  (33)  and 


f*(k,K,K)  =  ln(  —  \  -in  +  ttG  *  (42) 

V  |k| /  3 


SOLUTION  OF  THE  POTENTIAL  FUNCTION  WITH  THE 
SLENDER-BODY  ASSUMPTION 

Under  the  slender-body  assumption  the  length  of  the  ship  is  far  larger  than  the 
beam  or  the  draft.  We  separate  the  fluid  domain  into  two  regions:  the  outer  region 
where  (y,z)  is  of  the  order  of  the  length,  and  the  inner  region  where  (y,z)  is  of  the 
order  of  the  beam  or  the  draft.  In  the  outer  region,  the  three-dimensional  Laplace 
equation  is  solved  with  the  free-surface  condition.  Equation  (27),  and  the  radiation 
condition.  The  inner  solution  is  governed  by  the  two-dimensional  Laplace  equation 
with  the  free-surface  condition.  Equation  (28),  and  the  body  condition.  Equation 
(19),  on  the  ship  hull.  Then,  the  inner  and  outer  solutions  are  matched  in  the 
overlap  domain  to  determine  the  slender-body  solution. 
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Outer  Problem 


The  outer  solution  can  be  constructed  from  a  suitable  distribution  of  source 
strength  along  the  longitudinal  x-axis.  If  we  denote  the  source  strength  with  q^ (x)  , 
then  the  outer  solution  is  expressed  in  the  form 


=  J  Qj  (O  G(x-f,,y,z)  df 


for  j  =  3  and  5 


(A3) 


Here,  we  consider  the  solutions  for  heave  and  pitch  only.  The  Fourier  transform  of 
Equation  (43)  is  derived  from  the  convolution  theorem  in  the  form 


¥\*  =  qj*  G* (y , z ; k, k) 


(44) 


where  G*  is  defined  by  Equation  (39).  The  inner  approximation  of  Equation  (44),  for 
small  values  of  the  coordinates  (y,z),  can  be  constructed  by  substitution  of  Equation 
(41) 


ip  * 

j 


2D 


(43) 


Inner  Problem 

The  fundamental  solution  of  the  inner  problem  is  that  of  the  strip  theory. 
However,  the  matching  requirement  with  the  outer  solution  will  differ  from  the 
condition  of  outgoing  radiated  waves  satisfied  by  the  strip  theory.  Because  the 
outer  solution  includes  a  longitudinal  function  of  x  which  depends  upon  the  three- 
dimensional  shape  of  the  ship  hull,  we  solve  the  inner  solution  in  the  following 
f  orm^ 


J 


+  C (x) ) 


(46) 


where 


(s) 


Cj (x)  is  a  function 
is  the  strip  theory 


to  be  determined  by  matching  with  the  outer  solution  and 
solution  which  can  be  expressed  as  follows 
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(47) 


y  (s) 
j 


On  the  body  surface,  (J) .  and  cj>,  satisfy  the  following  conditions 


<j> .  =  -  iwn . 

jn  j 


(48) 


and 


(49) 


In  Equation  (46)  <t>  is  the  conjugate  of  <}>  .  The  solution  of  for  heave  motion  can 
be  expressed  by 

or 

A  _  _  „  ,  Y“\  T cos  2m 9  ,  K  cos  (2m-l)0  /rn 

*3  ~  a3G2D  \  I  2m  +  2^1  - 2^1 -  (5° 

m=l  r  r 

where  is  defined  by  Equation  (29),  is  the  two-dimensional  source  strength  at 
the  origin,  and  the  terms  under  the  summation  are  the  higher  order  multipoles  which 

form  the  wave-free  potentials.  Equation  (50)  was  first  introduced  by  Ursell."^  A 

/\ 

similar  solution  for  ^  can  be  given  by 


<t>3  °3G2d 


\  ^  [  cos  2m0 

am  I  2m 

m=l  r 


K  cos  (2m-l)0 


2m- 1 


2m- 1 


(51) 


For  pitch  motion  we  can  have  similar  expressions  to  Equations  (50)  and  (51)  by 

^  *  /N  /\ 

substituting  n  a  ,  o„,  and  a  for  -xa_,  -xa  ,  -xa_,  and  -xot  ,  respectively.  From 
JmJ  m  i  m  5  m 

Equation  (50)  the  outer  approximation  of  the  potential  <J>^  is 


=  °j  C2D(y’z) 


(52) 
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A  similar  expression  can  be  given  to  the  potential  <j>j  •' 


!>.  =  a.  G2D(y,Z) 


By  substitution  of  Equations  (52)  and  (53)  into  Equation  (46),  the  outer  approxi¬ 
mation  of  'f .  is 
1 


=  °j  G2D  +  °j  G2D  +  Cj(x)  (ajG2D+°jG2D) 


“  °j  G2D  +  °j  G2D  +  Cj(x)  (aj+Cj)  G2D  +  Cj(x)  (G2D-G2D)  °j 


From  Equation  (34),  the  last  term  of  Equation  (54)  is 


—  _  .  ,  Kz-iK  y  Kz+iK  y L  .  Kz  v  0. 

G2D  “  G2D  =  1^e  +e  )  =  2ie  cos  Ky  =  2i 


By  substitution  of  the  above  expression  into  Equation  (54),  the  outer  approximation 
is  finally  given  by 


V .  =  [o.+a.+C. (x)  (a.+o.)]  G„n  +  2i  C . (x)  a. 
J  JJJ  JJ  2D  j  j 


The  Fourier  transform  of  Equation  (55)  becomes 


V  .*  =  [o.+a.+C. (x)  (a.+a.)J*  G2D  +  2i[C.(x)o.]*  (56) 

J  111  J  J  31 


MATCHING 

The  inner  and  outer  solutions  are  matched  in  a  suitable  overlap  domain  to 

determine  the  unknown  source  strength,  q.,  of  the  outer  solution  and  the  coefficients  • 

C.  in  the  inner  solution.  From  Equations  (45)  and  (56)  we  have  the  following  1 

J  i 

relation  ! 
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(57) 


V  C2D  "  v  f*«*  =  [0j+Oj+Cj<°j-K,j)I*  G2D  +  21<Cjaj)* 


Equating  the  factors  of  G2^  gives  a  relation  for  the  source  strength 


q.*  =  [o ,+o ,+C . (o ,+o . ) ] *  for  j  =  3  and  5  (58) 

3  J  3  J  J  3 


Equating  the  remaining  terms  in  Equation  (57)  gives 


f  *q .  *  =  -2ni  [C.O. ]* 
3  3  3 


(59) 


The  inverse  Fourier  transforms  of  Equations  (58)  and  (59)  can  be  expressed  as 


q3 


[o ,+o  ,+C .  (o  ,+o . )  ] 
3  3  3  3  3 


2  tt  i  [C.o.  ] 
3  3 


f  (x-F,)dF, 


(60) 


(61) 


The  kernel  f (x)  is  the  inverse  Fourier  transform  of  Equation  (42).  Elimination  of 
C.  gives  an  integral  equation  for  the  outer  source  strength 


Vx> 


W 

2tto  . 

3  J 


I 


q.  (O  f(x-r,)dt;  =  (Oj+o.) 


(62) 


INNER  SOLUTION 

By  substitution  of  Equation  (61)  for  (x)  into  Equation  (46),  the  inner  so¬ 
lution  is  given  by 

<fi.  =  (<}>.+*.)  q.(0  f(x-F)dF  (63) 

33  2v~  3  3  J  3 

3  L 
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Equation  (63)  is  the  unified  inner  solution  which  is  valid  for  all  wave  numbers  K. 
The  details  of  the  derivation  and  error  estimations  are  given  in  Reference  3. 

In  Equations  (62)  and  (63),  when  j  is  5  for  pitch  motion,  the  two-dimensional 

'N  /X 

source  strengths,  and  a,.,  can  be  expressed  as  a  function  of  and  o^,  respec¬ 
tively-  The  two-dimensional  potential  for  pitch  motion  is 


*5  =  -  x  <5*3 


(64) 


*5  “  ^  *3  -  *  *3 


(65) 


From  these  equations,  0^  and  a ^  are  given  by 


n5  =  -  x  o3 


(66) 


°5  x  ''3 


ilT 


(67) 


By  substitution  of  Equations  (66)  and  (67)  into  F.qnation  (b2),  the  integral  equation 
for  the  outer  pitch  source  strength  is  given  by 


q5(x)  -  i 


(c5^5) 


2na 


3  J 


J' 


(][(")  fix-"  )<F  =  x('  +>’  )  +  -‘ 


iU 


(68) 


The  inner  solution  for  pitch  motion  becomes 


^5  =  ^  |q5(r)  f(x-r)df 

2tt(7  •' 


(69) 


The  derivation  of  the  kernel  function,  f(x),  from  f*(k,K,ic)  will  be  given  later. 
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DIFFRACTION  POTENTIAL 

With  a  similar  method  to  solve  the  potentials  for  j  =  3  and  5,  the  diffrac¬ 
tion  potential  can  be  solved.  In  the  outer  region  we  solve  the  three-dimensional 
Laplace  equation  and  in  the  inner  region  we  apply  strip  theory.  Through  matching, 
the  unknown  source  strength  and  coefficient  are  obtained.  In  the  application  of 
strip  theory,  however,  there  are  two  different  approaches.  One  is  to  solve  the  two- 
dimensional  Laplace  equation  and  the  other  is  to  solve  the  two-dimensional  wave 

equation  or  the  Helmholtz  equation.  The  first  approach  has  been  applied  by  Salvesen 
12  3  13 

et  al.  and  others.  Newman,  Troesch,  and  others  have  solved  the  Helmholtz 
equation  for  computing  the  diffraction  forces.  In  this  method  there  is  a  singu¬ 
larity  in  the  solution  when  the  angle  of  the  incoming  wave  is  0  or  180  degrees.  To 

3 

avoid  this  singularity,  Newman  has  introduced  the  long-wavelength  solution.  How¬ 
ever,  there  is  difficulty  in  solving  the  equation  for  short  waves. 

In  this  study  we  shall  solve  the  two-dimensional  Laplace  equation  to  compute 
the  diffraction  force  as  the  solution  of  the  strip  theory.  The  inner  solution  for 
the  diffraction  potential  can  be  given  in  the  form  (see  Reference  3) 

V  =  V^s)  +  C7(x)  (4-  +$  )  (70) 

/  /  /  s  s 


where  C^(x)  =  a  function  to  be  determined  by  matching  with  the  outer  solution 


yj  (s) 


=  the  diffraction  potential  solved  by  the  strip  theory 

(s) 


;  =  the  symmetric  function  of 

(  s  ) 

Tile  diffraction  potential,  *P  -j  ,  satisfies  the  two-dimensional  Laplace  equation  with 
the  boundary  condition  given  in  Equation  (18).  If  we  express  the  potential  of  the 
incoming  wave  as 


'P  =  -  — exp  (K  z+i  K  x  cos  P-i  K  y  sin  b) 
o  .'0  o  o  O' 

o 


the.  boundary 


cond it  ion 


of  P 


(s) 

7 


is  given  by 


(71) 


P 


(s) 
7  n 


K  (n  -i  n_  sin  ?)  P 
o  3  2  o 


(72) 
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r 


where  A  =  amplitude  of  the  Incoming  wave 

B  «  heading  angle  of  the  incoming  wave:  B  =  180  deg  in  a  head 
sea  and  8  *  0  deg  in  a  stern  sea 


K  =  oi7g 
o  o 


(») 


If  we  express  the  outer  approximation  of  Py  as 


=  n  C,  +  0 

7  7  2D  a  2D 


(73) 


then  the  outer  approximation  of  the  inner  solution,  Py,  which  is  similar  to  Equation 
(54)  ,  is  given  by 


^7  °7  G2D  +  C7  ^  G2D  +  C7(x^  ^G2D'G2D^  °7 


(74) 


The  second  term  in  Equation  (73)  is  asymmetric.  Because  the  ship  has  a  symmetric 
centarplane,  this  asymmetric  term  is  not  included  in  Equation  (74)  .  The  outer  so¬ 
lution  for  the  diffraction  potential  is  given  by  Equation  (45)  with  j  =  7.  Then, 
taking  the  Fourier  transform  of  Equation  (74)  and  matching  with  Equation  (45) ,  we 
can  derive  the  integral  equation  for  the  outer  source  strength  and  the  inner  so¬ 
lution  similar  to  Equations  (62)  and  (63)  as 


q?(x) 


I 


q 7  ( T )  f(x-OdF,  -  -  o. 


(75) 


V)  |Y)(s) 

7  7 


2tto. 


>/v 


f(x-r.)df 


(76) 


and 


C7  "  - 


-I  k 

2xa,  J 


(O  f(x-C)dr, 


(77) 
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HYDRODYNAMIC  FORCES 

The  pressure  in  the  fluid  is  given  by  Bernoulli's  equation 

P  =  -  P  (<t>t+gC+  \  Nl2^  (78) 

By  substituting  Equation  (3)  into  Equation  (78),  the  pressure  becomes 

P  =  p(iax|>2-V4yV<fr2)  e_iwt  -  Pgr.  (79) 


Then,  the  forces  and  moments  with  respect  to  the  origin  of  the  coordinate  system  are 
given  by 


F. 

J 


If  p  "j  ds 


for  j  =  3  and  5 


(80) 


where  S  is  the  submerged  portion  of  the  ship  hull,  and  n^  is  defined  in  Equations 
(20)  and  (21).  By  substitution  of  Equation  (79)  into  Equation  (80),  the  hydrodynamic 
forces  can  be  expressed  by 


F. 

J 


[  i£u<J)2-V<}>1*7i})0]  n.  dS 


-iwt 

e 


(81) 


The  second  term  of  Equation  (81)  can  be  transformed  by  means  of  a  theorem  due  to 
Ogilvie  and  Tuck  (Appendix  A  of  Reference  8) 


(V<j^-V<f>2)  n^  dS  = 


dS 


(82) 


In  Equation  (82),  the  line- integral  term  along  the  intersection  of  the 
the  plane  z  =  0  is  ignored  as  a  higher  order  term.  By  substitution  of 
(82)  and  (17)  into  Equation  (80),  the  hydrodynamic  forces  are  given  as 


ship  hull  with 
Equations 
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(83) 


-  - p  JJ  v  v.ZT5 


k=l 


dS  e 


-icot 


ADDED  MASS  AND  DAMPING 

We  write  the  integral  of  the  term  under  the  summation  in  Equation  (83)  as 


k  'k  6 


-160t 


(84) 


k=l 


where 


jk 


-  P 


If 


(  icon .  +m  .  ) 
J  3 


[<J'k+^k+Ck(x)(<})k+4>k)1 


dS 


(85) 


The  terms  inside  the  bracket  represent  the  unsteady  potential  f  which  results  from 

^  t-  Vi 

oscillation.  Equation  (85)  denotes  the  hydrodynamic  force  and  moment  in  the  j 
direction  per  unit  oscillatory  displacement  in  the  ktJl  mode. 

The  added  mass  and  damping  are  defined  by 


t 


jk 


i  B..  60 
jk 


(86) 


tor  j  =  j  ana  o  ana 


where  A.,  is  added  mass  and  B . ,  is  damping.  Both  A.,  and  u., 

jk.  J  k  jkjk 

k  =  3  and  5  are  given  in  Table  1.  In  Table  1,  is  used  the  following  notation: 


03  =  0R  +  fcfrj 

*3  =  K  +  **1 

C3  "  C3R  +  iC3I 
C5  "  C5R  +  iC5I 
n5  =  -xn3  ’  ^5  =  -x<^3 


m5  =  -(Un3+xm3);  iw  ^ 


=  777  ‘t’-i  “  x*. 
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The  added  mass  and  damping  in  Table  1  are  those  of  the  bare  hull  only.  The  added 
mass  and  damping  due  to  viscosity  and  stabilizing  fins  are  given  in  Table  2  and 
Table  3  of  Reference  1. 
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TABLE  1  (Continued) 


=  -  pJJ^R  dS  -  ^{[Vl  dS  "  Sjj^l  dS  "  n3^I  dS 

”  wjj  xm3<l>I  dS  '  P||  m3*R  dS  ~  2p  JJ  “^rSr  dS  ~  ¥  jj  n3RC3R 


^ff“3*RC3I  dS 


_r  S||  *2nA  ds  -  ff'H^i  ds  +  511 xn3  ds  ’  5If “3*1  ds 


^  2  JJ  3  yR  ^3  JJ  3  1 


+  ^2  JJ^R  dS  +  ^(f^R  dS  '  SJJx2n3^R  dS  "  f^jJ^R  dS 


+  ¥  [f  Xn3\C5R  dS  +  H 

(i >  JJ  OJ 


C5R  dS 


B5S  =  pjj  x2n3<PR  dS  +  ¥||  n3^R  dS  +  It||  Xn3*R  dS  +  ^2  JJ  “S^R  dS 


5|  x2m3^I  dS  +  ujj  X¥l  dS  +  P||  dS  +  2P||  ^VrSr 


+  ¥ffXn3W  ds  +  S*  ff  *S+RC5I  ds 


EXCITING  FORCES 

The  excitation  forces  are  given  by  Equation  (83)  as 


Fj  =  -p  ||  (iconj+mj)[  </>.,]  dS 


The  integral  containing  in  Equation  (87)  is  the  Froude-Krylov  force.  If  we 
substitute  <p  for  <f> ^  in  Equation  (82),  we  have 


■jj“j  ”0 da  *  -|/"Tr  ”j ds 


Here,  we  have  neglected  the  higher-order  terras.  With  Equation  (88),  the  Froude- 
Krylov  force  is  given  by 


f.  =  -P 


JJ  "J  (l 


ioH-U  ]  'fi  dS  e  1<l3t 

dX  }  O 


x)  ff  n.  * 

°  jj  j  ° 


where  ^  is  given  in  Equation  (71).  The  heave  excitation  force  is  given  by 


rr  KozT2 

f^  =  2pgA  JJ  1-e  |2  cos  (K^x-cut)  cos  K^y  dxdy 


and  the  pitch  moment  is 


rr  K^z  Z2 

f,_  =  2pgA  JJ  -1+e  0  ^  x  cos  (K^x-wt)  cos  K^y  dxdy 


+  2pgA 


—  fr e^°z (k  z-i)  2  sin  cos  K2y  dxdy 

K2  JJ  °Z_  Z1 

o 


where  K^  =  K^  cos  P,  sin  P,  and  the  integral  limits  z^  and  z^  are  given  in 

Figure  2. 
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Figure  2  -  Limits  of  Integration 

By  substitution  of  Equation  (70)  into  Equation  (87),  the  diffraction  force  is  given 
by 


h 

J 


(  icon ..  +nu  ) 


«P^s)+C7  (^s+4>s)  1  dS  e  1Wt 


(92) 


We  first  substitute  Equations  (48)  and  (49)  for  n.  and  m.,  and  apply  Green's  second 

(s)  J  J 


theorem  to  the  term  containing  *P 


The  diffraction  force  can  be  written 


h. 

1 


11  it  <VVds-°  11  (icon^-hiK  )  ((Jig+iJjg)  dS 


-io)t 


(93) 


(s) 

where  is  given  by  Equation  (77)  and  $  is  the  symmetric  part  of  -j  that  sat¬ 
isfies  the  two-dimensional  Laplace  equation  and  the  boundary  condition.  Equation 
(72).  By  substitution  of  Equation  (71)  for  ,  the  diffraction  heave  force  is 
expressed  by 


h^e 


iwt 


/  r  IKjX  r  KQz  \ 

li2poiriAje  dx  J  e  °  n^  cos  K2y-n2  sin  ®  *  sin  K2y)  (^3-03)  d(, 

c 

2p  C^(x)  dx  J  (ium^+m^)  Re[<f>g]  df.  (94) 

c 
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and  the  diffraction  pitch  moment  is 


,  icot  .  _  . 

h.e  =  -  i2pw  A 
5  o 


j"  (^x+  e"  ^  dx  J  e  °  cos  K^y-n^  sin  3  sin  ^y^  <f>3  d£ 


+  2p 


J  C-,  (x)  dx  J 


(iajn.x+Un-+m_x)  Re T 4>  ]  dS, 
3  3  3  s 


(95) 


Equations  (90),  (91),  (94),  and  (95)  represent  the  heave  force  and  pitch  moment 
of  the  bare  hull  of  the  SWATH  ship.  The  exciting  forces  and  moments  due  to  viscosity 
and  the  stabilizing  fins  are  given  in  Table  2  and  Table  3  of  Reference  1. 

EQUATIONS  OF  MOTION 

If  we  let  =  £3e~iwt,  a,.  =  £5e_1Wt,  the  total  exciting  forces  F^,  and  tht 
total  pitching  moment  F^,  then,  the  equations  of  motion  can  be  expressed 


(M+A.33)  0.3  +  B33CX3  +  C33oc3  +  A35a5  +  B35a5  +  C35ct5  =  F3 

( 1  ^ A 55)  “5  +  B55°l5  +  C55W5  +  A53^3  +  B53^3  +  C53a3  =  ^5 


(96) 


where  M  is  the  mass  of  the  ship  and  I  the  mass  moment  of  inertia  about  the  y-axis. 
The  hydrostatic  coefficients  are  given  by 


C33  Pg  Awp 


C_  _  =  C.,.  =  -  pg [ M  -x  • A  ] 
35  53  wp  g  wp 


(97) 


J55 


f  pg [ IT  +A  (x 
L  wp  g 


-xc)2]  + 


(KG-KB) 
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where  A  =  waterplane  area 
wp  r 

M  =  moment 
wp 

x  =  longitudinal  center  of  gravity 
8 

=  longitudinal  center  of  flotation 

1^  =  moment  of  inertia  of  A  with  respect  to  the  longitudinal 

center  of  flotation  ^ 

KG  =  distance  between  the  keel  and  vertical  center  of  gravity 
KB  =  distance  between  the  keel  and  vertical  center  of  buoyancy 
A  =  displacement  of  the  ship 

With  substitution  of  =  F33  e  l0Jt  and  F<.  =  F^  e  ia)t,  Equation  (96)  can  be 
expressed  as  complex  algebraic  equations: 


0)  (M+A33)  +  C33  -  iwB33  u  A35  +  C35  -  iojB35  I  C3 

-  0)2  A53  +  C53  -  iu>B53  -  a32(H-A55)  +  C55  -  LuB^  I  £5 


RELATION  OF  TWO-DIMENSIONAL  SOURCE  STRENGTHS  BETWEEN  THE  FRANK  CLOSE-FIT 
METHOD  AND  MULTIPOLE  EXPANS T CM  METHOD 

The  potential  function  given  in  Equation  (5 (j)  is  the  solution  obtained  by  the 
multipole  expansion  method.  If  we  apply  the  Frank  close-fit  method  to  solve  the 
two-dimensional  potential,  we  should  have  some  relationship  between  the  two  methods. 
That  means  we  should  have  a  correspondence  between  the  sources  located  at  the  origin 
in  the  multipole  expansion  method  and  those  around  the  ship's  contour  in  the  Frank 
close-fit  method. 

Because  the  radiation  condition  states  that  far  away  from  the  disturbance,  the 
wave  is  outgoing,  we  should  have  the  same  potential  from  these  two  methods  at  a 
distance  far  from  the  body.  For  the  multipole  expansion  method,  the  potential  of 
the  outgoing  wave  is  given  by  Equation  (34) 


=  -io  e 
o 


Kz+iK I y  I 


where  a  is  the  source  -strength  at  the  origin, 
o 
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For  the  Frank  close-fit  method,  we  have  the  potential  as  follows 


c 


(100) 


where  a  is  the  complex  source  density  at  the  ship's  contour  and  G.  for  the  twin-hull 

9 

body  is  given  by  Wehausen  and  Laitone 


G2D  “  Re  2tt 


log  (y+iz-n-iO  -  log  (y+iz-n+iO  +  2PV 


1 


=-ik(y+iz-ri+i^) 

K-k 


dk 


-  i  Re  [e-iK(^+iz-ri+i*’>  ] 


+  Re 


2tt 


log  (y+iz+n-i?)  -  log  (y+iz+n+iC)  +  2PV 


3-ik(yHiz+ri+iC) 

K^k 


dk 


-  i  Re  [e-iK(y+iz+T1+i^)  ] 


(101) 


where  (y,z)  is  the  point  where  the  potential  is  sought  and  (n,C)  is  a  source  point 
(Figure  3) 
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when  | y |  is  large,  the  logarithmic  terms  cancel  out.  Then,  Equation  (101)  can  be 
expressed  in  a  manner  similar  to  Equation  (30) 


2D 


oc 


gk(z+0  cog  k(  n) 

- iUk - -  dk  + 


OC 


k(z+0  -  ,  . 

® - c°s  him l  dk  d02) 

K-k 


Equation  (102)  can  be  further  transformed  with  the  exponential  integrals  as  shown  in 
Equation  (31).  As  | y  j  becomes  large,  these  exponential  integrals  become  small. 
Finally,  as  jy|  becomes  large.  Equation  (102)  can  be  expressed  as 


=  aK(-z+o  r ^ iK < y_n '  +  eiKiy+ri 


C2D  =  6 


(103) 


Then,  the  potential  function  is  given  by 


=  |  o[-ieK(z+rJ  eiKiy-n']  de 


1 


+  I  n[-ieK(z+r')  e1K|y+n!]  d« 


.  .  Kz  iKv 
!  i  e  e 


cos  Kr,  df 


(104) 


Equating  Equation  (99)  with  Equation  (104),  we  have  the  following  relation  between 

•  and 
o 


a 

o 


2 


KC 

e  cos 


Kn  df. 


(105) 


27 


This  equation  simply  states  that  the  corresponding  source  at  the  origin  is  obtained 
by  integrating  the  distributed  source  along  the  contour  with  the  weighting  function 
of  2e^  cos  Kn. 


DERIVATION  OF  THE  KERNEL  FUNCTION 

The  inverse  Fourier  transform  of  Equation  (42),  we  can  express  f  (x)  as'* 


f (x)  =  £n(2K)  5(x)  +  -  +  tt  G,,(x,0,0)  -  irr  5(x) 

2 1  x  I  3 


(106) 


where  t5(x)  is  Dirac's  delta  function  and  is  defined  by  1  where  x  =  0,  and  0  when 
x  #  0.  From  Equations  (C3  through  C4)  and  (35),  G^  (x,0,0)  is 


(x,  0, 0)  =  ~  I  e  U 
4:t“  J 


J  77  ?  2 

(k“+!!“)  [(k~+?z)  -k] 


r  rv  1/2  l/9  “ 

1  I  (K  ~~4m  ~u  cos  [1)  exp(-ixu  cos  0) 

2lT  o  o  u  -  (K1/2+m1/2u  cos  0) 


dO  (107) 


where  m  =  U"/g. 

The  poles  of  Equation  (107)  are 


7  \  1/7  1/7 

\  1  -  2(mK) _ cos  0  J  [l-4(mK)  cos  0 

i,  I  2  m  cos^  0 


The  integral  path  of  Equation  (107)  is  determined  by  the  value  of  cos  "  as  follows: 
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1/2 

CASE  1  -  4(mK)  cos  0  <  1 

The  integral  limit  in  0  is  between  6^  and  n/2,  where  6^  is  given  by 


0  =  cos 

T 


-1 


MinK)1/2 


=  COS 


-i/JJ 

4t  I 


The  integral  path  of  u  is  given  in  Figure  4. 


1  L1 
u-> 


►  U 


1/2 


Figure  4  -  Integral  Path  when  4(mK)  cos 


CASE  2  -  4(mK)1/2  cos  0  >  1 

The  poles  in  Equation  (108)  become 


Ul\  1/2  1/2 

1  cos  '■  *  i[4(mK) _ cos  0  -1] 

,  2 
2  m  cos  (• 


The  integral  limit  in  ''  is  between  0  and  . 

CASE  1  -  cos  0  <  0 

With  the  change  of  the  variable  0  =  n  -  C- ,  the  integral  with  respect,  to 
becomes 


I  :  ^  -  I 

tt/2  o 


1/2  1/2  —  2  — 

(K  -m  “u  cos  0)  exp(ixu  cos  0)  ^ 


lvl/2  1/2 

u  -  (K  -m  u  cos  ■ ) 


(109) 


(110) 


dll) 


The  poles  of  Equation  (111)  are 


A  1/9  _  1/2  _ 

4  '  1  +  2 (mK) 1/ ~  cos  9  ±  [ 1+4 (mK)  '  cos  6] 

2  — 

2  m  cos  6 


(112) 


The  integral  path  of  u  is  given  in  Figure  5. 


i 


Figure  5  -  Integral  Path  when  cos  0  <  1 
By  substitution  of  Equation  (111)  into  Equation  (107),  C, ^  (x,0,0)  is 


C3(x,0,0) 


1/2  1  /  0 

(K  '  A  cos  0)  exp(-ixu  cos  9) 

u  -  (K^2+m1/2u  cos  6) 


+ 


1/9  1/0 

(K  A  cos  •')  exp  ( ixu  cos  ”) 

du - - 9 - 

,,,1/2  1/2 

u  -  vk  -m  u  cos  ) 


(113) 


With  the  following  expansion 


(K 


,1/2  i/: 


0) 


K  +  u  cos 


I/O 

[mu  cos  0  *  2 (mK)  “ ] 


Equation  (113)  can  be  rewritten 
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G  i  f  d6  du  +  r  de  I  du  -K-  exP(-ix-u  cos  0)  , 

2tt  \  J  J  {  I  ,„l/2  1/2  ‘ 

\o  o  0  /  u-(K  +m  u  cos  0) 


M  de  I 


K  exp(ixu  cos  6) 


2tt  J  j  1/2  1/2  _  ‘ 

o  u-(K  -m  u  cos  0) 


J  00 


r4I  da  / 


du  +  d0  du 


p  o  0  L . 

<  T  i 


u  cos  P[mu  cos  0+2 (mK)  ]  exp(-lxu  cos  ') 
u-(K1//24m1,/2u  cos  0) 


£  J 


d0  I  du 


o  Ln 


u  cos  0 [mu  cos  0-2 (mK)  ‘ 

,VU2  1/2  ‘ 

u-(K  ~m  u  cos  0) 


exp(ixu  cos  0) 


g3  ~  Dx  g4 


(114) 


where  and  are  given  by 


g3  2 


<*  T  n/2  . 

J  dO  j  du  +  j  dfi 


K  exp(-ixu  cos  0) 


L}  /  u-(K1/2+m1/2u  cos  0) 


TT  /  2  - 

?  f  * 


K  exp(ixu  cos  0) 


j  i  fvU2  in  ■ 

o  I,  u-(K  -m  u  cos  0) 


-I?(f  de  1  du+{  deJ  d“)- 

0  x  1 1  / 


1/2 

-i[mu  cos  Q+2(mK)  ]exp(-ixu  cos  8) 
u-(K1/2-Hn1'/2u  cos  0) 


(116) 


*/2  , 

£  J  d9{  d“‘ 


1/2 

[mu  cos  0-2 (mK)  }exp(lxu  cos  0) 
u-(K1//2-m1//2u  cos  0) 


With  the  change  of  the  integral  paths  and  the  change  of  the  variables,  we  can 
derive  the  equations  suitable  for  the  numerical  computation  from  Equations  (115)  and 
(116).  Following  the  derivation  of  Joosen,^  we  can  express  and  as  follows 
(see  Appendix  D) : 


■  -  8(1> 
’n  ®n 


+  g 


(2) 


(3) 


for  n  =  3  and  4  (117) 


where 


(1)  _ 


Hi 


2(2)  it 


'  ,  2  1/2  1/2  1/2  1/2 

Rn(v)<[  [(m^+m0)  +m^ ]  +i(sign  n^)  [  (m^+n^)  -m^] 


-v  x 


2  2 
(m^+m9) 


1/2 


—  dv 


(118) 


(2) 


2ik 

IT 


4J 


K  or  -i 


+2  (mK)  1  ^ 


L\ 

f[k1k+2(mK)1/2]  -4k2k2} 


1/2 


exp 


kk 

-i  — —  x)  dk  for  x  <  0  (119) 
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(2) 


2ik„ 


a 

I 


K  or  -i 


^2^  1/2 


71  J  i/o  4  ?  2  ly^ 

1  {[k^+ZCmK)1^]  -4^} 


eXP  (-1  2^X)dk 


21k, 


a 

I 


K  or  i 


k4k  1/2 

-  -2(mK)I/Z 


IT  J  1  /2  4  ?  ? 

1  {[k4k-2(mK)  '  ]  -4k^k  / 


1/2  £XP  (*  2m-  X)  dk 


2ik„ 


1 

I 


K  or  i 


k3k  1/2 

L~~  -2(mK)i/2 


*o  C[k3k-2(mK)  -4kZkZ} 


^  N^(k)  exp  |  -  —  [1-2  (mK) 


(3)  =  I  f  In 

*  J 


k^k 

L/j  exp  2^—  x  )  dk  for  x 


1/Z]  xjexp  (-  ^  6kx) 


0  (120) 


1/2 


[k(l-k)  (Sk+1)  (6k+2) ] 


exp  j-(5(l-k)  (Sk+1) ] 


1/2 


2m 


(121) 


In  Equations  (119)  and  (120),  the  first  numerator  is  for  n  =  3  and  the  second  for 
n  =  4.  All  other  notations  used  in  Equations  (118)  through  (121)  are  as  follows: 


R3(v)  =  K 


1/2 

R  (v)  =  -i  (mK)  -  mv  sign  x 


=  (K-mvZ)  -  (4  mK-1)  vZ 


(122) 
(cont  .  ) 


1/2  2 

=  4  (mK)  v  (K-mv  )  sign  x 
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k  .  1  -  2(mK)1/2  -  [l-4(mK)xr"l 


.1/2 


1/2 


=  2  (mK) 


1/2 


k?  =  1  -  2(mK)1/2  +  1 1—^* (niK^ "l’/  *"] 


,1/2 


1/2 


=  2  (mK) 


1/2 


=  1  +  2  (mK)1^2  -  [ 1+4 (mK) 2 ^ 2 ) 


1/2 


k  =  l  +  2  (mK)1/2  + 
4 


1/2 


1/2 


6  =  4(mK)1/2  -  1 


for  (mK)^2  <  0.25 


for  (mK)^  >  0.25 


for  (mK)^2  <  0.25 


for  (rnK)^  >  0.25 


(122) 


N3(k)  -  K 


N  (k)  =  -  --  [6 (1-k)  (5k+l)]1/2  Sign  X  -  ±  [6k+l+2(mK)A,i] 


1/2. 


K  = 


*  in  special  cases,  when  K  »  0,  G3  becomes  (see  Appendix  D) 

Ho  (m)+  (2+  Sign  X>  Yo(t) 


1 

C3  =  4  Ox 


(123) 
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and  when  m  -  0 , 


(124) 


G,  =  7  [-H  (K | x  j )  -  Y  (K | x I )  +  21  J  <K | x  J  > ] 

J  h  o  o  o 

where  H  is  the  Struve  function,  and  J  and  Y  are  Bessel  functions  of  the  first  and 
o  o  o 

second  kind. 


RESULTS  AND  DISCUSSION 

To  test  the  numerical  results,  we  have  selected  three  models  of  SWATH  ships: 
SWATH  hA,  SWATH  6C,  and  SWATH  6D.  Experimental  tests  for  motion  have  been  carried 
out  for  all  three  models.  All  models  have  the  same  lower  hulls  and  the  same  distance 
between  the  two  hulls,  but  differ  in  the  shape  of  their  struts.  SWATH  6A  has  a 
single  strut  on  each  side  while  SWATH  ships  6C  and  6D  have  twin  struts.  The  thick¬ 
ness  of  the  struts  differs  for  all  three  models.  The  principal  dimensions  of  these 
models  are  given  in  Table  2.  Because  our  primary  interest  is  the  motion  at  low  fre¬ 
quencies  of  encounter,  the  computations  have  been  carried  out  when  the  forward  speed 
is  20  knots. 


TABLE  2  -  PRINCIPAL  DIMENSIONS  OF  SWATH  SHIPS* 


Dimension 

SWATH  6A 

SWATH  6C 

SWATH  6D 

Displacement,  long  tons 

2579 

2602 

2815 

Length  at  waterline,  m 

52.5 

60.2 

68.0 

Length  of  main  hull,  m 

73.2 

73.2 

73.2 

Beam  of  each  hull  at  waterline,  m 

2.2 

2.6 

2.2 

Hull  spacing  between  centerline,  m 

22.9 

22.9 

22.9 

Draft  at  midship,  tn 

8.1 

8.1 

8.1 

Maximum  diameter  of  main  hull,  m 

4.6 

4 . 6 

4 . 6 

Longitudinal  center  of  gravity 
aft  of  main  hull  nose,  m 

35.5 

34.7 

36.1 

Vertical  center  of  gravity,  m 

10.4 

10.4 

9.  1 

Longitudinal  CM,  m 

6.1 

13.7 

26.4 

Radius  of  gyration  for  pitch,  m 

16.9 

17.7 

19.0 

Waterplane  area,  m^ 

193.9 

181.4 

253.9 

Length  of  strut,  m 

52.4 

25.8  (per 
strut) 

25.8  (per 
strut ) 

Strut  gap,  m 

0 

8.6 

16.4 

Maximum  strut  thickness,  m 

2.2 

2.6 

3.  1 

*Dlmensions  are  full-scale. 
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Figures  6,  7,  and  8  show  the  added  mass  and  dampirg  coef f icients .  The  presently 

computer  results  are  compared  with  those  of  strip  theory.  The  results  of  strip 

theory  were  computed  by  the  computer  program  MOT35  (Reference  16)  which  is  based  on 

the  analytical  method  developed  by  Lee.^  The  two-d impnsional  potential  (strip 

theory)  is  solved  by  the  Frank  close-fit  method.  The  effects  of  fins  and  viscosity 

are  included  in  the  computation  of  added  mass,  damping  forces,  and  excitation  forces. 

This  computer  program  has  been  improved  by  adding  surge  effect  to  the  pitch  excita- 

2 

tion  moment  due  to  the  incoming  waves  and  by  correcting  for  viscous  damping. 

The  results  for  added  mass  by  the  present  method  are  almost  70  to  80  percent 
higher  than  those  from  strip  theory,  while  the  damping  coefficients  are  slightly 
less  than  those  from  strip  theory.  In  the  present  method  the  results  of  the  unified 
slender-body  theory  are  added  to  those  of  the  strip  theory  as  a  correction  term 
(Equation  (63)).  For  the  added  mass  coefficients,  these  correction  terms  are  too 
large.  The  large  discrepancy  may  result  from  the  different  methods  used  to  solve 
the  strip  theory.  In  the  derivation  of  the  unified  slender-body  theory,  Newman 
solved  the  strip  theory  with  the  multipole  expansion  method  (Equation  (50)).  In 
Equation  (50)  a  source  located  at  the  origin  and  is  ushd  to  solve  the  three- 

dimensional  source  strength  (Equation  (62)),  from  which  the  correction  factor  Ch  is 
computed  as  in  Equation  (61) .  The  only  part  of  the  source  located  at  the  origin  is 
o^.  The  other  part  of  the  source  is  from  Equation  (50),  which  is  included  in  the 

computation  of  the  two-dimensional  added  mass  and  damping  coefficients,  but  is  not 
included  in  the  three-dimensional  source  strength. 

In  the  present  approach  for  solving  the  two-dimensional  problem,  the  Frank 

close-fit  method  (Equation  (100))  was  applied.  In  Equation  (100)  a  is  distributed 

on  the  contour  of  the  section,  and  the  corresponding  coerce  a  ,  which  is  located 

at  the  origin,  is  computed  by  the  relation  in  Enuf-ion  (105).  To  solve  the  three- 

dimensional  source  strength,  was  replaced  with  o  ,  but  Oq  includes  the  effect  of 

o _  and  a  .  It  is  difficult  to  decompose  a  into  o_  and  a  in  the  Frank  close-fit 
3  m  o  3  m 

method.  If  cr^  is  not  computed  accurately,  the  three-dimensional  source  q^  i:i 
Equation  (62)  cannot  be  solved  correctly,  even  though  the  kernal  function  f (x)  has 
been  evaluated  properly.  Therefore,  the  difference  between  and  Oq  could  possibly 
cause  the  large  discrepancy  in  added  mass  coef f icients. 
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STRIP  THEORY 


igure  f>  -  Added  Mass  and  Damping  Coefficients  of  SWATH  6A  in  Following 

Seas  at  20  Knots 


Figure  8  -  Added  Mass  and  Damping  Coefficients  of  SWATH  6D  in  Following 

Seas  at  20  Knots 


1/2 

When  (mK)  =  o)U/g  =  0.25,  there  is  a  discontinuity  in  the  computations.  This 

1/2 

singular  behavior  can  be  explained  by  the  Green  function.  If  (mK)  approaches 
0.25,  Equations  (119)  and  (120)  become  logarithmic  functions  as  shown  in  Reference 
15. 

Figures  9  and  10  show  the  results  of  excitation  forces  and  moments  for  SWATH 
6A  and  6D,  respectively.  The  heave  forces  for  both  configurations  are  slightly 
larger  than  those  from  strip  theory.  The  excitation  moments  are  generally  larger 
also.  For  SWATH  6D,  the  results  of  strip  theory  for  the  heave  force  show  better 
agreement  with  experimental  results  than  those  from  the  unified  slender-body  theory. 
The  pitch  moments  of  SWATH  6D  are  scattered  between  the  results  of  the  strip  theory 
and  the  experiment. 

The  motion  results  for  SWATH  6A  are  plotted  in  Figure  11.  Heave  amplitude 
shows  a  similar  discrepancy  with  experiment  as  with  strip  theory.  When  the  en¬ 
counter  frequency  becomes  very  small,  the  peak  of  the  heave  amplitude  disappears; 
this  peak  is  shifted  to  the  higher  encounter  frequencies.  The  plus  points  show  the 
results  of  the  "mixed  method,"  in  which  the  added  mass  and  damping  coefficients  are 
computed  by  the  unified  slender-body  theory  and  the  excitation  forces  are  computed 
by  the  strip  theory.  When  A/L  is  less  than  2.5,  heave  amplitude,  computed  by  this 
method,  agrees  better  with  experimental  results  than  with  the  results  of  either  the 
unified  slender-body  theory  or  the  strip  theory.  However,  when  A/L  is  larger  than 
2.5,  the  results  of  the  mixed  method  underpredict  the  experimental  results  and 
those  of  both  theories. 

Pitch  amplitude  shows  the  same  tendency  as  heave.  The  peak  is  shifted  to 
higher  A/L  values.  Pitch  amplitude,  computed  by  the  mixed  method,  agrees  well  with 
experiment  for  all  A/L  values. 

The  motion  results  for  SWATH  6C  are  given  in  Figure  12.  The  heave  and  pitch 
amplitudes  become  very  large  as  the  encounter  frequency  becomes  large.  The  results 
of  the  strip  theory  are  better  than  those  of  the  unified  slender-body  theory  and 
the  mixed  method.  From  the  results  of  the  mixed  method,  we  can  conclude  that  the 
added  mass  and  damping  coefficients  computed  by  the  unified  slender-body  theory  in 
Figure  7  are  too  large. 
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Figure  10  -  Exciting  Force  and  Moment  of  SWATH  6D  in  Following 

Seas  at  20  Knots 


X 

L 


I  i.iin.  il  -  Motion  ut  SWA  IH  f»A  in  Following  Soas  at  J  0  Knots 
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Figure  13  shows  the  computed  motions  of  SWATH  6D.  The  heave  amplitudes  pre¬ 
dicted  by  tfie  slender-body  theory  are  much  larger  than  those  predicted  by  the  strip 
theory.  Compared  with  experiment,  the  results  of  the  mixed  method  are  best  except 
for  high  */L  values. 

The  pitch  amplitudes  computed  by  the  unified  slender-body  theory  and  by  the 

mixed  method  are  too  large  and  do  not  show  good  comparison  with  the  results  of  strip 

theory  or  the  experiment.  In  contrast  to  the  results  for  SWATH  6A,  the  mixed  method 

does  not  compute  the  pitch  amplitude  properly.  This  indicates  that  the  added  mass 

and  damping  coefficients  in  Figure  8  by  the  slender-body  method  may  be  it  significant 

error.  This  error  might  he  caused  by  solving  the  strip  theory  with  the  Frank  close- 

fit  method  md  by  replacing  0 _  with  o  as  mentioned  in  the  discussion  of  the  added 

3  o 

mass  and  damping  coefficients. 

SUMMARY  AND  CONCLUSIONS 

3 

In  this  report  the  unified  slender-body  theory  developed  by  Newman  is  applied 
to  predict  the  motion  of  SWATH  siiips  in  following  seas.  Only  for  a  limited  range  of 
•/I,  values  is  there  an  improvement  for  heave  motion.  For  pitch  motion,  except  small 
encounter  frequencies,  the  results  are  worse  than  those  of  strip  theory.  When  the 
encounter  frequency  is  large,  the  pitch  notion  becomes  extremely  large.  The  reason 
for  this  discrepancy  seems  to  lie  in  the  fact  that  the  strip  theory  is  solved  by  the 
Frank  close-fit  method,  ami  not  by  tire  multipole  expansion.  From  the  present  study, 
the  following  conclusions  may  he  drawn: 

t.  in  the  unified  slender-body  approach,  the  largest  contribution  to  the 
Lvdrodvnani  ic  roof  f  io  ients  results  from  the  solution  of  the  strip  theory  itself. 

rr<.<t  computation  of  the  source  strength  at  the  origin  is  necessary  to  compute  the 
t  hree-d imens ional  source  strength  exactly.  Therefore,  for  the  outer  approximation 
of  strip  tlieorv,  the  multipolo  expansion  method  should  be  applied  instead  of  the 
Frank  close-fit  method. 

A  more  careful  analysis  should  be  made  in  computing  the  excitation  forces 
bv  the  strip  theory.  Application  of  the  Helmholtz  equation  is  mathematically  more 
correct  .  However,  there  is  a  singularity  in  this  solution  when  the  heading  angle 
of  the  incoming  waves  is  zero  or  180  degrees.  The  two-dimensional  Laplace  equation 
.(imputes  the  heave  excitation  force  correctly,  but  the  pitch  excitation  moment  is 
com  [Hi  t  od  i  nciirriT  t !  y . 
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1.  According  to  the  unified  slender-body  theory  it  should  be  valid  for  all 
t requeue  ies.  However,  for  large  frequencies,  the  motion  results  become  unacceptably 
large.  This  may  be  caused  by  the  use  of  a  different  method  for  the  solution  of  the 
strip  theory. 

a.  Because  of  the  above  mentioned  discrepancies,  until  further  research 
mentioned  in  items  1  and  2  are  carried  out,  the  application  of  the  unified  slender- 
body  theorv  to  improve  the  prediction  of  motion  of  SWATH  ships  is  not  recommended . 
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APPENDIX  A 

DERIVATION  OF  G2D(y>z)  WITH  EXPONENTIAL  FUNCTION 


Equation  (30)  can  be  rewritten  as 


G2D(y,z) 


-i-  f 

2  it  J 


,kz(eiky+e~lky) 


k-K 


dk 


(A.  1) 


For  the  first  term  of  Equation  (A.l),  we  change  the  variable  as 


w  =  -i(k-K)  (y-iz) 


(A. 2) 


and 


dk 


-dw 

i(y-iz) 


(A.  3) 


The  contour  of  the  integral  path  will  be  different  depending  upon  y  <  0.  When  y  >  0, 
we  take  the  integral  path  in  the  following  figure 


IP.  substituting  Equations  (A. 2)  and  (A. 3)  into  Equation  (A.l),  the  first  term 

!>«.••■  i  >!'!!•  . 


r  kz+iky 

J  k-K 


I 


iK(y-iz)  -w 
e  e 


dw 


iK(y-iz) 


(A. 4) 
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If  we  consider  the  following  integral  along  the  closed  path  in  the  w-plane 


J Vd“-J+  J  +J- 

^7  x  r  v 


2tt  i 


then,  the  integral  over  f  =  Re^  vanishes  as  R  +  and  Equation  (A. 4)  becomes 


^  dk  __  eiK(y-iz) 
k-K 


iK(y-iz) 


2Tii-  |  - dw 

w 


=  elK(y  1Z)  [27Ti+E1(Kz+iKy)  ; 


(A. 5) 


When  y  <  0,  the  integral  path  is  given  in  the  following  figure 


k-PLANE 


and  the  first  term  of  Equation  (A.l)  becomes 


oo  ro-|-^oo 

C  kz+iky  f 

|  j 

~  -tr/ 


iK(y-iz)  -w 
e  e 


dw 


(A. 6) 


iK(y-iz) 


If  we  consider  that  the  following  integral  along  the  closed  path  is  in  the  w-plane. 
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v  d»  •  I + 1  *  S  ■ 
^  ^  r  - 


then  the  integral  over  T  =  Re^  vanishes  as  R  +  and  Equation  (A. 6)  becomes 


I‘ 


£2*  dt .  _ 


k-K 


J 

J 


iK(y-iz) 


dw 


eiK(y-iz)  Ei(Kz+iKy) 


(A. 7) 


For  the  second  term  of  Equation  (A.l),  we  apply  the  following  transformation  of  the 
variable 


w  =  i(k-K)  (y+iz) 


(A. 8) 


and 


dk 


dw 


i(y+iz) 


(A. 9) 


Cith  the  same  procedure  for  derivation  of  Equations  (A. 5)  and  (A. 7),  the  second  term 

>f  Equut  ion  (A.l)  is  given  ns 


r  kz-iky  _ 

J  — -£■  -- —  dk  =  e  Z  1  y  E^(Kz-iKy)  for  y  >  0  (A. 10) 

o 


=  eKz  [  2iri+E^  (Kz-iKy)  ]  for  y  <  0  (A. 11) 


Bv  substitution  of  Equations  (A. 5),  (A. 7),  (A. 10),  and  (A. 11)  into  Equation  (A.l), 

< ,  is  ;■  i von  bv 

.’It 
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APPENDIX  B 


DERIVATION  OF  G,D(y,z)  FOR  SMALL  ARGUMENT 

We  expand  the  exponential  function  and  the  exponential  integral  in  Equation  (31) 
for  small  Kr  where  y  =  r  sin  0  and  z  =  -r  cos  0 

eKz+iKy  =  1  +  Kz  +  iKy(l+Kz)  (B.l) 

^Kz+iK  |  y  |  =  j  +  Kz  +  iKjyJd+Kz)  (B .  2) 

E1(Kz+iKy)  =  E  (Kre1  *  ) 

=  -y  -  iJn(Kr)  — i (tt— G )  +  Kz  +  iKy  (B.3) 

By  substitution  of  Equations  (B.l),  (B.2),  and  (B.3)  into  Equation  (31),  is 

given  as 

G2n(y,z)  =  -  ^  Re  { [  (1+Kz+iKy  (1+Kz)  ]  [-y-£n(Kr)-i(iT-0)+Kz+iKy]  } 

-  i  [ (l+Kz)+iK j v i (1+Kz) ) 

=  —  (1+Kz)  [  y+?n  (Kr ) -Kz+Kv‘l- i”  ]  (B.4) 
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APPENDIX  C 

FOURIER  TRANSFORM  OF  THE  THREE-DIMENSIONAL  GREEN  FUNCTION 

j 

Equation  (36)  can  be  expressed  as  i 


G(x,y,z)  =  - 


4  it 


/ 


-ikx 


dk 


J 


°°  (k2+£2) 


1/2 


2  2  1/2 
,  ,  k  exp[z(k  +8.  )  -iy£]  I 

1/2  ‘ 

2  2 

(kZ+IT)  -  K 


dSt  (C.l) 


where 

K  =  (C .  2) 

g 

The  first  double  integral  of  Equation  (C.l)  is  known.  Therefore,  we  rewrite  Equation 
(C.l)  as 


G(x,y,z) 


1 

2ttR 


G3(x,y,z) 


(C.3) 


2  ?  2  2 

where  R  =  x“  +  y  +  z  and  G^  is  given  as 


C3(x,y,z) 


k  I  Ah  J 


2  2  1/2 

K  exp[z(k  +8.  )  — iyg. ] 

2  2  1/2  22  1/2 
-»  (k  +8.  )  [(k  +8.  )  -k] 


d£ 


}  (C.4) 


If  we  assume  y  and  z  become  very  small.  Equation  (C.3)  is  expressed  as 


G(x,0,0) 


-  G,(x,0,0) 


2  2 

27t(x  +r  ) 


1/2  3 


(C.  5) 


.  2  2  2 
where  r  =  v  +  z  . 


iriLUL.ui.vC 


LLjuJC-iiiOT  FI 
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The  Fourier  transform  of  Equation  (C.5)  is 

G*(k,K)  =  -  ^  Ko(|kjr)  -  G3*(k,K) 


(C.6) 


and  from  the  definition  given  in  Equation  (38),  G^*  is  given 


G3*(k,K) 


1  I  k  dS. 


J  2  2  1/2  22  1/2 
(k  +S,  )  f  (k' +1  )  -k] 


(C.  7) 


In  Equation  (C.5),  K  (x)  is  the  modified  Bessel  function.  For  small  r,  K  (Iklr)  can 

o  o 

be  expanded 


K0([k[r)  =  -  y 


(C.8) 


where  y  =  0.577...  is  Euler's  constant. 

By  substitution  of  Equation  (C.8)  into  Equation  (C.6),  for  small  y  and  z,  G*  is 


C*(k,K) 


g3* 


=  G 


2D 


f *(k,K,<) 


(C .  9) 


where  is  Equation  (33),  G3*  is  Equation  (C.7),  and  f*  is  given  by 


f*(k,K,K)  =  J,n  -  -in  +  ttG,* 

|ki  3 


(C .10) 
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where  and  are  given  in  Equation  (110).  The  first  integral  in  Equation  (D.l) 
is  the  same  as  the  first  in  Equation  (115).  The  second  integral  along  the  contour 
c^  becomes  zero 


1 


(  )  dw  =  lim 

R-KC 


-tt/2 

J 


i0 


KjZZizl |g_  x-  -C°—  9)  Re19  i  d9  -  0 
Re10  -  (  )2 


The  last  integral  of  Equation  (D.l)  becomes 


r  c 

J  (  )  dw  =  j* 


K  exp(-iivx  cos  8) 


n  1 


iK  exp(-vx  cos  8)  dv 


(  )  l  i/2  1/2  2J  1/2  1/2  ^ 

iv  -  [ (K)  +(m)  iv  cos  0]  o  iv+[(K)  -i(m)  v  cos  0] 


Therefore,  the  first  integral  of  Equation  (115)  is  given  by 


w  — 

C  t  i  q \  exp(-iu0x  cos  9) 

f  - — J  du  ..  2* - - - 2 - m 

o  u  -  I (K) ^^+(m) cos  0]  [4(mK)^^  cos  6  -1] 


j 


iK  exp(-vx  cos  8) 


dv 


(D.2) 


o  iv  +[  (K)  1/,2-i(m)  ^^2v  cos 


b.  Second  Integral  of  Equation  (115) 

Evaluate  the  following  integral  in  the  complex  plane  of  w 


K  exp(-iwx  cos  9) 


w  -  [  (K)^2+(m)  ^2w  cos  6] 


dw 


■S  ■!  ■!■ 


2iri  (residue) 
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The  residue  is 


K  exp(-iu2x  cos  0)  -K  exp(-iu?x  cos  9) 


2  a/  x  ,  /0  1/2 

°(u2  V  [l-4(mK)1/2  cos  0] 


The  integral  along  c^  becomes  zero  as  shown  in  (a)  and  the  integral  along  becomes 


J  (  )  dw  =  J  - 


K  exp(-j.ivx  cos  6) 


-00  iv-  [  (K)  '*'^2+i(m)  cos  0] 


idv 


-J: 


iK  exp(-ivx  cos  6) 


dv 


o  iv+  [  (K)  ^2-i(m)  ^‘'v  cos  0] 


Therefore,  the  second  integral  of  Equation  (115)  is  given  by 


) 


(  )  du  =-  2?rKi 


K  exp(-iu0x  cos  ft) 


[l-4(mK)^2  cos  0 


or 

J 


iK  exp  (-vx  cos  0) _  dv  (D>3) 


^2  J  j/ 2  1/2  2 

|  o  iv+[ (K)  ~-i(m)  v  cos  6] 
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c.  Third  Integral  of  Equation  (115) 

Evaluate  the  following  integral  in  the  complex  plane  of  w 


r 

O 


_ K  exp(iwx  cos  6) 

..,,,1/2  ,  ,1/2  ‘ 

w-[(K)  -(m)  w  cos  @] 


dw 


2Tii  (residue) 


K  exp(iu^x  cos  9)  K  exp(iu^x  cos  0) 

2  +  2 
-m  cos  "(u.-u.)  -m  cos'-  0(u_-u.) 

4  3  3  4 


-K  exp(iu^x  cos  r')  K  exp(iu^x  cos  0) 


l/o 

1+4  (mK)  *"  cos 


1/2 


1/7 

1+4 (mK)  "  cos  0 


1/2 


The  integral  along  c9  becomes  zero  as  shown  in  previous  cases  and  the  integral  along 
c^  becomes 


1  ‘  ’  d"  ’  1 


_ K  exp(iivx  os  0) 

1/2  1/2 

iv-[(K)  -(m)  iv  cos  0. 


idv 


o 


iK  exp (-vx  cos  ‘ ■)  dv 
1/2  l/o 

iv-[(K)  -i(ra)  ^v  cos  9] 


o 
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r  K  exp(iu_x  cos  ■') 

J  7/T-  Tri  d  ' 

*o  [l+4(mK)  '  cos  6  ] 


2.  x  <  0 

Wo  can  follow  the  same  method  to  transform  the  integral  with  respect  to  u.  But, 
in  this  case,  we  take  different  integral  paths:  for  the  first  and  second  integrals, 
we  eviluate  the  complex  integrals  given  in  Case  1  in  the  first  quadrant,  and  for  the 
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third  integral,  we  evaluate  the  complex  integral  in  the  fourth  quadrant.  The  ex¬ 
pression  for  g^  similar  to  Equation  (D.5)  is  given  by 


ti/2 

■  *  /  '  J 


_ K  exp(vx  cos  0)  dv 


l/9  i/-> 

o  iv-[(K)  “+i(m)  “v  cos  0] 


7 ' 

£  I  J"  J  . 


iK  exp (vx  cos  0)  dv 


l /->  [  /  •> 

o  iv+[(K)  ~+i(m)  “v  cos  0] 


J 


K  exp(-iu.x  cos  0) 


o  [  4  (mK)  ^^cos  0-1] 


1/2 


dB 


,  r 


K  exp(-iUjX  cos  0) 


O  [1-4  (mK)  '  ^ “  cos  0  1 


i/: 


dd 


(D.6) 


For  g,  ,  the  same  procedure  can  be  applied.  But,  the  numerator,  K,  should  be 

1/2  1/9 

replaced  by  -i[mu  cos  0+2 (mK)  “|  or  i[mu  cos  0-2 (mK)  " ]  with  appropriate  variables 
and  poles. 

I f  we  let 

A3  =  K 


B3  =  K 


.sl/2, 


(D.7) 


A,  =  -i[mu  cos  0+2 (mK)  “ ] 

*  i/-» 

B,  =  i [mu  cos  0-2(mK)  ~ ) 

we  can  express  g^  for  n  =  3  and  4  as  follows 


r  f  A  (  -  in » o) 
do  j  - 


A  (  -  i  u  )  exp(-ux  cos  0) 


du 


1/1 


o  iu+[(K)  "-i(m)  "u  cos  ('  ] 


j  -n  / 


B  (iu,0)  oxp(-ux  cos  0)  du 
1/’  l/’ 

iu-[(K)  ^  —  i  (in)  "u  cos  0] 


1 


A  (i^.O)  exp(-iu0x  cos  0) 


1/2 

[ 4 (mK)  cos  0-1] 


1/2 


d0 


(D.8) 

(  r  on  r  .  'I 
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tt/2 


+  1 

71 


I 


An(u2>P)  exp (- iu^x  cos  0) 


"■/ 2 


[ 1-4 (mK) ^ ^ cos  0 1 


1/2 


dO 


•  /•  *-» 
«  I 


B  (u,,u)  exp(iu,x  cos  ") 


1/2 


1/2 


o  [l+4(mK)  “  cos  ! 


tt/2 


i  f  _J1_3 

71  J 


B^(u.},0)  exp/iu^x  cos  0) 


,1/1 


2/2 


dO 


[l+4(mK)  “  cos  0] 


for  x  >  0 


(D.8) 


tt/2 


■'£  I  d7 


A  (iu,0)  exp(ux  cos  0) 


du 


1/9  1/9 

o  iu-[(K)  ‘'+i(m)  "u  cos  0] 


tt/2 


£  I  d°  I 


B  (-iu,0)  exp(ux  cos  0) 


du 


1/9  1/9 

o  iu+[(K)  “+i(m)  “u  cos  0] 


-if  ^ 

IT  0 


A^(u,,0)  exp(-iu^x  cos  0) 


1/2 


dO  + 


[ 4 (mK) 


1/2 


cos  0-1) 


i7 


tt/2 


A  (u.,0)  exp(-iu  x  cos  0) 

- - - jy2—  do  (D .  9) 

[ 1-4 (mK) ^  cos  0] 


for  x  <  0 


The  first  two  double- integral  terms  in  liquations  (D.8)  and  (D.9)  can  be  reduced 
with  substitution  of  u  cos  0  =  v  (Reference  17)  to: 

1.  n  =  3  and  x  >  0  f or  =  B.^  =  K 
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i 


I 


=  (sign  m^)  r  sin  0 


cos  6  = 


(m2-hn2) 


1/2 


sin  0  = 


lm2  I 


,  2  2. 

(m  +m2) 


1/2 


the  numerator  becomes 


1/2  1  22  1/2  1/2  2  2  1/2  1/2 
(ra^+in^)  =  ~Yj2  f(mi+m2^  +miJ  +i(sign  m2>  [ (m^+n^)  -m^] 


Finally,  can  be  expressed 


(1) _ L 


2(2)  1/2tt 


/ 


2  2  1/2 
^ml+m2^  +ml 


+i(sign  m.) 


r  2  2  1/2  11/2 1 

(m^+mp  -m^ 


-vx 


(m2+m2) 


1/2 


dv 


(D. 12) 


2 .  n  =  3  and  x  <  0 

Following  the  same  method  as  that  of  Case  1,  we  can  derive 


(1) 


2  (2)  1/2tt 


on 

f 

r  r  2  2  1/2  i 

1/2 

2  2  1/2 

1/2 

J 

(m^+m2)  +m 

+i(sign  m9) 

(m^+m9)  -m^ 

O 

VX 


(m2+m2) 


1/2 


dv 


(D. 13) 


where 


4(mK)^2  v  (K-mv2) 


(D. 14) 
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With  Equations  (D.12)  and  (D.14),  can  be  written  for  all  values  of  x  as 

follows 


(1) 


2  (2)  1/2tt 


1/2 


+m. 


1/2 


+i(sign  m2) 


,  2  2. 

(m^+n^) 


1/2 


1/2 


-v  x 


(m^-Hn^) 


1/2 


dv 


(D.15) 


where 


m2  =  4(mK)'*'^  v  (K-mv^)  (sign  x) 


(D. 16) 


When  n  =  4,  by  following  the  same  procedure,  the  double-integral  terms  in 
Equations  (D.8)  and  (D.9)  can  be  expressed  for  all  values  of  x 


.(1) 


CO 

t175;  J 


2(2)x/^ 


R4(v) 


2  2^2 

(m1+m2>  +m1 


1/2 


+i(sign  m2) 


(m^+m2) 


1/2 


1/2 


-V  X 


2  2 
(u^-Hi^) 


1/2 


dv  (D. 17) 


where 


R4(v)  =  -i  2(mK) 


1/. 


-  mv  sign 


(D.18) 


The  rest  of  the  terms  in  Equations  (D.8)  and  (D.9)  can  be  further  reduced  to 
the  forms  that  are  more  useful  for  numerical  computations. 

For  x  <  0,  let  the  last  term  of  Equation  (D.9)  be 


tt  /  2 

r  An(U 

(2)  _  i 

n  " 71 i  fl 


A^u^O)  exp(-iu^x  cos  0) 


[1-4 (mK) cos  q] 


1/2 


d0 


(D . 19) 
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where 


A3 (u  1 , 0 )  =  K  and  A^u^.O)  =  -itmu^ 


1/2 

cos  0+2 (mK)  ] 


With  the  change  of  the  variable 


2  mu^  cos  0  =  k^k 


(D. 20) 


where  is  given  in  Equation  (108)  and 


k  =  1  -  2(mK)1/2  -[l-4(mK)1/2]  for  (mK)1/2  <  0.25 


=  2 (mK) 


1/2 


1  /? 

for  (mK)i/4  >  0.25  (D.21) 


(2) 


is  transformed  as 


i^i  _  c 

TT  J  - 


K  or  -i 


^1^  1/2 
+2(mK)i/Z 


1/2  2  2 
[k  k+2(i«K)  '  ]  -4k^k 


j  1/2 


exp 


.  V  \ 

1  2m  */ 


dk 


(D . 22) 


In  the  numerator  of  Equation  (D.22),  the  first  term  is  for  n  =  3  and  the  second  for 
n  =  4 . 

For  x  >  0,  to  the  fourth,  fifth,  and  sixth  terms  of  Equation  (D.8),  we  sub¬ 
stitute  the  following  variables 


2  111U2  cos  0  =  k£k 


2  mu,  cos  8  =  k,k 
4  4 


2  muj  cos  8  =  k^k 


(D. 23) 


(D.24) 


(D. 25) 
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The  third  term  of  Equation  (D.9)  for  x  <  0  becomes 


i |T  Vv 


8)  exp(-iu^x  cos  8) 


o  f4(mK)  cos  8-1] 


-1 


2(mK) 


1-2 (mK) 


K  or  |-|  (4mK-k2)  -i  j  +2  (mK)  1/2J  j  dk 
1/2  ^[2(mK)1/2-k][k+2(mK)1/2]  [k+2(mK) 1,2  ]  -1  j 


ikx  .  r.  ,  1/2  , 2,  x 

exp  {-  -r-—~  +  [  4  (mK)  -k  ]  -r- 


(D. 31) 


The  third  term  of  Equation  (D.8)  for  x  >  0  can  be  transformed  to  the  form  similar 
to  Equation  (D.31).  Because  of  the  lengthy  derivation,  we  omit  the  intermediate 
steps. 

With  the  change  of  the  variable  in  Equation  (D.31) 


t  m  k  2(mK)1/2-l 

K  1/2  "**  1/2 
4(mK)  -1  4 (mK) 1  -1 


(D . 32) 


can  be  finally  expressed  as  follows  for  all  values  of  x 
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’n  it  J 
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(D . 33) 
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where 


6  =  4(mK)1/2  -  1 

N3  =  K  (D. 34) 

N4  =  -  \  [6(l-k)(6k+l)]1/2  (sign  x)  -  |  [ 6k+l+2 (mK) 1/2 ] 

Then,  and  g^  given  in  Equations  (115)  and  (116)  can  be  expressed  with  n  =  3 

and  4 


(1) 

!n  =  gn 


+  g 


(2) 


(3) 


(D. 35) 


where  g  is  Equation  (D.15)  for  n  =  3  and  Equation  (D.17)  for  n  =  4;  ^  is 

n  (3)  n 

Equation  (D.22)  when  x  <  0  and  Equation  (D.29)  when  x  >  0;  and  g^  '  is  Equation 

(D.33).  For  the  details  of  these  derivations,  see  Reference  15.  In  this  reference, 

there  are  some  errors  in  signs. 

As  in  special  cases,  when  m  =  0  or  K  =  0,  we  can  integrate  these  equations 
explicitly . 

1.  Steady  Forward  Motion,  K  =  0 
Because  of  K  in  the  numerator. 


g 


(1) 

3 


0 


From  Equations  (D.10),  (D.16),  and  (D.18), 

2/  2  2. , 
=  v  (m  v  +1) 


m^  =  0 

R^(v)  =  -  mv  sign  x 
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By  substitution  of  these  relations  into  Equation  (D.17),  becomes 
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1/2  e  1  1 
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(D. 36) 


In  Equation  (D.36),  H  (x)  is  the  Struve  function  and  Y  (x)  is  the  Bessel  function  of 

o  '  o 

the  second  kind. 

(2) 

When  K  =  0,  k^  =  k^  =  0  and  k^  =  k^  =  2.  Therefore,  for  x  <  0,  =0. 

(2) 

For  x  >  0,  from  Equation  (D.29),  is  given  by 
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(D- 37) 


71 


T 


(3) 

The  integral  limit  of  in  Equation  (D.31)  is  transformed  from  the  variable  0 

that  is  between  0  and  0  =  cos  ^  {l/[4(mK)2//2]} .  If  1/ [4  (mK)  ^2  ]  is  larger  than  1, 

T  (3) 

O  becomes  zero.  Therefore,  when  K  =  0,  g^  becomes  zero.  By  adding  Equations 
(D.36)  and  (D.37),  is  given  by 


g4 


-  \  H 
4  o 


{2+  sign  x]  Yq 


Finally,  from  Equation  (114),  G3(x,0,0)  is  given 
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(D.38) 


(D .39) 


2.  Pure  Oscillation,  m  =  0 

Because  R,  ,  Equation  (D.18),  is  zero,  gf^  becomes  zero.  When  m  =  0,  k.  =  k~  = 

’  (2)  ^  ^  J 

0  and  k?  =  k^  =  0.  Therefore,  g^  ,  Equation  (D.22),  is  zero  when  x  <  0.  When 

x  >  0,  as  the  complex  argument  of  the  exponential  function  becomes  infinitive  in  the 

(2) 

first  and  second  terms  in  Equation  (D.29),  g,  becomes  zero.  The  upper  limit  of 

(3)  4 

the  integral  for  g^  in  Equation  (D.33)  origined  from  9^  which  becomes  zero  in  this 
case.  Therefore,  there  is  no  contribution  from  g^  to  when  m  =  0. 

With  R3  =  K,  m3  =  K2  +  v2,  and  m,  =  0,  g^  is 
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is  evaluated 


With  A3  =  K  and 
from  Equation  ( 


Because  =  K, 


For  all  values 


Because  of  0 


u  =  K  from  Equation  (115),  as  m  becomes  zero, 
.9)  for  x  <  0 
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for  x  <  0  (D .41) 


and  u0  and  u^  disappear. 


for  x  >  0  is  given  from  Equation 


(D.8) 
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for  x  >  0  (D.42) 


of  x,  is  given  by 
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exp(iK | x 
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=  |  [-HQ(K|xi )+iJ0(Kix| ) ] 


(D.43) 


0,  g3  is  zero .  Finally,  g3  is  g 


iven 


G  =  e  =t  [-H  (K , x ' ) -Y  ( K , x | )+2iJ  (K | x | ) ] 
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1.  Page  20,  Table  1:  Equation  for  A^  at  bottom  of  table,  replace  fourth  term 


*  £  JJWr  ds 


jjn3*R  dS 


2.  Page  21,  Table  1  (continued):  Equation  for  A,.,.  at  middle  of  table,  replace 
eighth  term 

'  S£//*2n3*R  dS 


IF" 


3*R  C5I  dS 


3.  Page  35,  Table  2:  Value  given  for  "Waterplane  area"  for  "SWATH  6D"  as 
"253.9,"  fourth  line  from  bottom  of  table,  should  be  211.2. 

4.  Page  51,  top  of  page:  Insert  integral  sign  in  front  of  quantity  on  left- 
hand  side  of  equation  so  that 


becomes 


r  -W 

Jv 


